Scattering and absorption of ultracold atoms by nanotubes 
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We study cold atoms scattering from static and vibrating nanotubes. We find the scattering po- 
tential calculations for dielectric tubes can be greatly simplified with a modified pairwise summation 
approach. Quantum reflection of atoms from a nanotube is shown to become important if a Bose- 
Einstein condensate is used. Non-trivial absorption of atoms by nanotubes, as seen in experiments, 
is explained by inter-atomic interactions and quantum pressure. Natural nanotube vibrations do 
not significantly reduce condensate fractions, but lower frequency oscillations can dramatically heat 
the cloud. 

PACS numbers: 34.20.-b, 34.50.Cx, 61.48.De, 67.85.-d 
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Interfacing ultracold atom clouds with micro- and 
nanoscaled objects enables exciting possibilities for cre- 
ating hybrid quantum systems which combine the long 
quantum coherence times of ultracold atoms with the 
flexibility of other molecular, optical and solid state tech- 
nologies PHI]- There has been work on cold atoms 
with ions and superconductors 0-Q > proposals for hybrid 
nanosystems [7H9j and recent experiments in which the 
interaction between cold atoms and carbon nanotubes 
was observed and measured Tot 11|. 

However, there is much to learn about how atoms scat- 
ter from nanostructures, or are absorbed by them. This 
is because theories of dispersion forces only become sim- 
ple in very basic geometries where the finite size of the 
objects is neglected; this often cannot be done in the 
case of nanosystems. In addition, scattering of ultracold 
atoms is expected to take place in the crossover region be- 
tween the non-retarded and retarded zones of the poten- 
tial, which complicates the theory further. Recent exper- 
iments with ultracold atoms and carbon nanotubes raised 
further questions regarding the applicability of classical 
scattering theories 12 1, non-dynamical approaches 
the role of quantum reflection and the effect of vibrations 
in the nanotubes. Better understanding of the dispersion 
interactions and scattering in such systems is crucial to 
future hybrid experiments and will aid the design of nan- 
otechnologies where such forces exert a profound influ- 
ence on their mechanical properties. It may also provide 
insight into the manufacture of nanostructures since this 
is often done from the vapor phase. 

In this paper, we study how cold atoms are either 
elastically scattered from, or absorbed by, a nanotube. 
We begin by calculating the Casimir-Polder potentials 
for conducting and insulating tubes and then use these 
potentials to study the likelihood of elastic scattering by 
considering quantum reflection probabilities. We then 
perform dynamical simulations of a Bose-Einstein con- 
densate (BEC) interacting with a nanotube to determine 



the rate at which atoms are lost due to inelastic scat- 
tering. Finally, we consider how an oscillating nanotube 
can affect the condensate fraction of a cold atom cloud 
when the two are overlapped, again by performing time- 
dependent simulations. 

The interaction between an atom and a cylindrical, 
electrically neutral solid object is expected to be dom- 
inated by dispersion forces, specifically the Casimir- 
Polder (CP) potential. This subject has been consid- 
ered by several groups 13- H) and is complicated by the 



fact that the properties of the cylinder alter not only the 
strength coefficients, but also the functional form of the 
potential [l3^ . 

The potential in the case of an atom and an infinitely 
long tube can be calculated using the general Lifshitz 



formula 
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V CP (r) = /i fc S T^'£ ; 2 c^)Tr \ G ^(r,r,i^)\ (1) 



y 2 is the separation between the atom 



where r = yfx 2 
(in the x-y plane) and the tube axis z, T is the tem- 
perature of the tube, = 2TrksTl /h are the imaginary 
Matsubara frequencies with mode number I (the dash on 
the sum implies half weight for the first term), a(i^i) is 
the polarizability of the atom, and the trace is taken over 
G( s )(r, r, the scattering part of the Green's polariz- 
ability tensor for the tube. 

Equation ([1]) is appropriate for conducting and insu- 
lating tubes, but the potential in the insulating case can 
be approximately modeled by using the simpler pairwise 
summation approach, which is used for AFM tips and 
nanostructures [22] and is efficient for arbitrary geome- 
tries. The interaction energy can then be expressed as 



Vcp{t) = - I ^^du <hui: 



(2) 



V pws 



where r pws is the separation between the atom and a 
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FIG. 1: (a) CP potential energy curves for a nanotube. Con- 
ducting tube, Liftshitz model (gray curve), dielectric tube, 
pairwise summation model (blue curve) and the — Cs/r 5 po- 
tential (red curve). Inset: effective power law, k, for the three 
cases, color scheme as in the main plot. The dotted line marks 
the position of the tube surface. 



discrete piece of the tube, C(r, £) is the strength coeffi- 
cient, and the integral takes place over the volume of the 
tube, V. C may be calculated as a function of r and £ to 
provide approximate corrections due to retardation and 
temperature effects. We use [l^| 
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where r/ = (2^/c)r ptos with c the vacuum speed of light. 
/3(i£i) is the polarizability of a discrete piece of nanotube 
used for the pairwise summation. As a(i£) generally de- 
cays much more rapidly with frequency than /3(«£), we 
can use static tube polarizabilities to a good approxi- 
mation. Finally, we consider the simplest approxima- 
tion Vcp{r) — —Cs/r 5 which is a limit of the pairwise 
approach, assuming L 3> r 3> Rq m the non-retarded 
regime, where Rq is the radius of the tube. 

We plot the potential for 87 Rb atoms interacting with a 
nanotube at T = 300 K in Fig. [1] assuming R = 63.5 nm 
[Ill |. The 87 Rb polarizability is taken from [H]. We 
plot one curve (gray) using Eq. fl} and the Green's ten- 
sor from [2l|, which is for a good conductor. The blue 
curve is plotted using Eqs. ((2j) and ([3]) with a tube length 
L = 10 (im and j3 calculated assuming a dielectric con- 
stant e = 2.5. The red curve is plotted using the —C$/r 5 
model with = 6 x 10~ 65 Jm 5 , as obtained in recent 
experiments [lOj. We see that the potential becomes 
strong on ultracold atomic energy scales (~ 10~ 31 J) 
about 0.5 /im from the surface. The conducting tube 
generally delivers the strongest potential, as expected. 
The —Cs/r 5 potential is around three orders of magni- 
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FIG. 2: Quantum reflection probabilities, R, from the dif- 
ferent nanotube potential models. Conducting tube (gray 



curve), dielectric tube (blue curve) and the 
tial (red curve). 



-Cs/r poten- 



tude weaker. We calculated the potential for a dielec- 
tric cylinder using the Green's tensor approach [24[ but 
have not plotted it here since it is barely distinguishable 
from the pairwise summation curve, the difference being 
~ 5%. This is significant since Eq. ^ can typically be 
calculated orders of magnitude faster than Eq. ([T]) and 
is immediately extendable to arbitrary geometries. It is, 
however, crucial that retardation corrections are included 
using Eq. 

The inset of Fig. [1] shows the effective power laws for 
the curves, k — d(log Vcp)/d(\og r), plotted against r. 
The gray curve (conducting tube) goes to n = —3 at 
the surface and tends to k w -4 with increasing r. By 
contrast, the dielectric model (blue curve) is closer to 
k = —5, indicated by the red line. 

We now analyze quantum reflection of atoms from the 
tube. For simplicity, we restrict our analysis to the per- 
pendicular component of incident A:-vectors, k perp and 
reduce to problem to one dimension (ID). The quan- 
tum reflection probabilities can then be calculated with 
ID Schrodinger wavepackets in a harmonic trap (25| . 
The condition for strong quantum reflection is given by 
$(fc) = (l/k 2 )dk/dr - d\ dB /dr > 1 [26] where the local 
wave number, k — 2ir/\dB, depends on the de Broglie 
wavelength, A<js, for the atom's center-of-mass motion. 
This defines a quantum reflection zone, which typically 
lies between r ~ 60 nm and r ~ 2 /jm for cold atoms. 
Atoms which come closer than the inner boundary can- 
not reflect and are absorbed by the nanotube. We model 
this effect with an imaginary potential [27|. The reflec- 
tion probabilities, R, for the three different potentials 
are plotted in Fig. [2] for perpendicular incident velocities 
in the range v perp — 5 mm s which 

are common to ultracold atom clouds. We see that the 
conducting tube (gray curve) barely reflects 87 Rb atoms 
at the relevant ultracold velocities, which is expected 
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for stronger potentials [28|. The insulating tube reflects 
atoms better but only delivers probabilities above ~ 0.2 
when the velocities drop below ~ 2 mm s _1 . This corre- 
sponds to a thermal cloud with T ~ 40 nK. The — C5/V 5 
potential delivers noticeably higher reflection probabili- 
ties but they are still <C f for clouds with T > 1 00 nK, 
mean speed > 5rnms . We note that quantum reflec- 
tion is higher for the nanotubes than would be expected 
for planar surfaces since the potentials are steeper 
and better fulfill the quantum reflection condition. 

The results suggest that quantum reflection is likely 
to play an important role if a BEC is made to interact 
with a nanotube. A stationary BEC has no center-of- 
mass motion, but the curvature of the wavefunction can 
create a quantum pressure on an absorbing surface, which 
leads to atom loss, and hence, atom flow with velocities 
on the order of < lmms" 1 (2f| . We now investigate 
the interaction of a 87 Rb BEC with a nanotube. We use 
the simple — C5/V 5 potential and assume a cylindrical 
symmetry with the tube fully overlapped in the center 
of the BEC. We model the dynamics of the condensate 
using the Gross-Pitaevskii equation (CPE) 
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where a = 5.4 nm is the s-wave scattering length, V 2 is 
the Laplacian in cylindrical coordinates, and ^l{z,T,t~) is 
the wave function at time t, normalized such that \^\ 
is the number of atoms per unit volume. The external 
potential, Ve, is a combination of the —C§/r potential 
and the harmonic trap. As with the reflection simula- 
tions, we employ imaginary potentials to absorb atoms 
that do not reflect. We choose lj z = 2tt x 80radss _1 
and u r = 2ir x 37 rads s" 1 in the vertical and horizontal 
(radial) directions respectively, which are close to recent 
experiments (lQ| . The initial atom number is chosen by 
averaging the first two points of the experimental curve 
we wish to com par e with. Spin flip losses are expected 
to be negligible 21| . 

We plot the atom number against time in Fig. [3] for 
the simulation (blue curve) and the corresponding ex- 
perimental loss curve (black points). The loss curve is 
non-trivial showing a wide range of loss rates but the 
agreement with experiment is good. Atoms are accel- 
erated towards the tube and the slow initial loss rate 
increases as atoms are rapidly removed from the BEC. 
This generates collective oscillations in the trap, shown 
in the inset of Fig. |3l along with sounds waves and phase 
gradients that can produce interference patterns. The 
subsequent fast loss rate flattens into a plateau; at this 
point very few atoms are near the trap center due to the 
in-trap oscillations. Finally for t > 0.05 s we see a rel- 
atively slower loss rate; at this point most of the atoms 
have been absorbed so the repulsive inter-atomic inter- 
actions, and hence the quantum pressure are weaker. 
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FIG. 3: BEC losses due to a nanotube, atom number plotted 
against time. Black points are experimental data, blue curve 
is for the GPE simulation with interactions, magenta curve is 
for a non-interacting wavepacket. Inset show the cross-section 
of the BEC's density profile after 0.03 seconds. Gray shapes 
sketch the relative position of the tube and the chip surface. 



For comparison, we also show a simulation with no 
inter-atomic interactions [Fig [31 magenta curve]. All 
other parameters remain the same. We see a classic expo- 
nential decay, with a loss rate that is significantly slower 
than in the experiment. This reveals that the interactions 
are important, since they increase the quantum pressure 
of the wavefunction on the tube. We note that the very 
simple model for the CP potential is sufficient to describe 
the experimental data. 

Up until now we have assumed the nanotube may be 
treated as a static object. However, it is possible that 
the interaction between atoms and nanotubes is signifi- 
cantly affected by thermal or mechanical oscillations of 
the nanotube. This might produce deleterious effects but 
conversely, might prove useful in cooling nanoscaled ob- 
jects. We therefore end our investigation by considering 
the effects of a vibrating nanotube on a BEC. This time 
we consider a full three-dimensional system with the po- 
tential further simplified to a step function. The dynam- 
ics are again modeled by the Gross-Pitaevskii equation 
with 10 4 87 Rb atoms, trap frequencies as before. The 
simulation begins with the BEC ground state in the trap 
and the tube potential is ramped up over 1 ms which sim- 
ulates insertion of the tube. The tube is then oscillated 
sinusoidally for 0.02 s with an amplitude of A = 1.5 /Ltm, 
variable frequency, /. This may be achieved by driving 
the tube with a piezo element. Spatial and energy cut- 
offs absorb atoms that are excited to very high energies 
by the process; we typically lose 80-90% of the atoms. 
The tube potential is then suddenly turned off and the 
one-particle density matrix of the cloud is then evalu- 
ated after averaging over 0.04 seconds. This reveals the 
relative occupation of the ground state mode, / ci after 
Penrose- Onsager [29| . 

We plot f c against the vibration frequency of the tube 
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FIG. 4: Remaining condensate fraction against nanotube os- 
cillation frequency following a 0.02s overlap. Insets show 
cross-section cloud profiles at the end of the oscillation pro- 
cess. 



in Fig. [U The insertion and removal of the tube reduces 
f c to about 0.6. We see that at low frequencies this is 
largely maintained. However, at higher frequencies, the 
tube oscillations reduce f c to around 0.2. At even higher 
frequencies, f c returns to a relatively constant 0.6 (the 
complex dynamics lead to small fluctuations in the ex- 
tracted value) . 

This behavior is explained by the BEC's speed of 
sound, v s , and correlation time, t c = Ih/v s where Ijj is 
the condensate's healing length. Topological excitations 
and reduction of f c first set in when the speed of the tube 
is greater than the local speed of sound, A ■ 2irf > v s . 
This corresponds to / > 60 Hz in this case and the pop- 
ulation of excited modes increases with the frequency. 
However, if the tube oscillates with a frequency / 3> l/i c , 
the gas cannot respond to the oscillations and merely sees 
a static, smeared potential; the disruption is minimal, 
and f c is retained [Fig. [3] inset bottom right]. In our 
case, l/t c ~ 3000 Hz, assuming the peak atom density. 
Between these limits, the BEC is strongly disrupted [Fig. 
IH inset top left]. 

Nanotubes typically have a fundamental frequency 
^ 100 kHz. We therefore conclude that natural vibra- 
tions of nanotubes are unlikely to disrupt the BEC, in 
spite of the high energy of the oscillations. Furthermore, 
undriven thermal fluctuations for typical nanotubes have 
an amplitude < 5 nm [3fj| . This implies the tube would 
not move faster than the speed of sound in the BEC and 
very little disruption would occur. 

In conclusion, we have studied quantum reflection from 
nanotubes and shown it becomes significant for Bose- 
Einstein condensates. We saw that the nanotube poten- 
tial extracted from recent experiments was significantly 
weaker than that expected for metallic tubes or even 



a good insulator. This might be explained by surface 
roughness effects [25, 28]. Pairwise theories are only guar- 
anteed to become accurate in the limit of rarified media 
but we find a close match with the Lifshitz theory in 
the dielectric case; this could potentially make design of 
nanodevices much simpler. 

We then simulated BEC atom losses to a nanotube, 
finding highly non-trivial dynamics. A detailed descrip- 
tion of the wavepacket, including quantum pressure and 
inter-atomic interactions was necessary to describe the 
experimental data. However, only very simple theories 
of the potential were required. Interestingly, no other 
loss channel, other than Casimir-Polder physics, was re- 
quired to explain the results. Finally we considered how 
an oscillating nanotube excites a BEC. We found there is 
a resonance where the condensate fraction is significantly 
reduced at a frequency < l/t c . Very high and very low 
frequencies damage the condensate much less. 
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